Skin, as the outmost layer of human body, is frequently exposed to environmental stressors including pollutants and ultraviolet (UV), which could lead to skin disorders. Generally, skin response process to ultraviolet B (UVB) irradiation is a nonlinear dynamic process, with unknown underlying molecular mechanism of critical transition. Here, the landscape dynamic network biomarker (l-DNB) analysis of time series transcriptome data on 3D skin model was conducted to reveal the complicated process of skin response to UV irradiation at both molecular and network levels. The advanced l-DNB analysis approach showed that: (i) there was a tipping point before critical transition state during pigmentation process, validated by 3D skin model; (ii) 13 core DNB genes were identified to detect the tipping point as a network biomarker, supported by computational assessment; (iii) core DNB genes such as COL7A1 and CTNNB1 can effectively predict skin lightening, validated by independent human skin data. Overall, this study provides new insights for skin response to repetitive UVB irradiation, including dynamic pathway pattern, biphasic response, and DNBs for skin lightening change, and enables us to further understand the skin resilience process after external stress.
